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Abstract 

Directional detection of Galactic Dark Matter is a promising search strategy 
for discriminating WIMP events from background. Technical progress on 
gaseous detectors and read-outs has permitted the design and construction 
of competitive experiments. However, to take full advantage of this powerful 
detection method, one need to be able to extract information from an ob- 
served recoil map to identify a WIMP signal. We present a comprehensive 
formalism, using a map-based likelihood method allowing to recover the main 
incoming direction of the signal and its significance, thus proving its Galac- 
tic origin. This is a blind analysis intended to be used on any directional 
data. Constraints are deduced in the (a n ,m x ) plane and systematic studies 
are presented in order to show that, using this analysis tool, unambiguous 
Dark Matter detection can be achieved on a large range of exposures and 
background levels. 
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1. Introduction 



A substantial body of astrophysical evidence supports the existence of 
non-baryonic Dark Matter : locally, from the rotation curves of spiral galax- 
ies Pi, on cluster scale from the observation of galaxy cluster collisions j^j 
and on the largest scale from cosmo logical observations [3]. Candidates for 
this class of to-be-discovered particles (WIMP: Weakly Interacting Massive 
Particles) naturally arise from extensions of the standard model of particle 
physics, such as supersymmetry [IJ or extra-dimensions j^]. 
Like most spiral galaxies, the Milky Way is supposed to be immersed in a 
halo of WIMPs which outweighs the luminous component by at least an order 
of magnitude. Tremendous experimental efforts on a host of techniques have 
been made in the field of direct search of WIMP Galactic Dark Matter. With 
really small expected event rates (O(10~ 5 — 1) day _1 kg _1 ), the challenge is 
to distinguish a genuine WIMP signal from background events. Two main 
Dark Matter search strategies may be identified : the detector may be de- 
signed to reach an extremely high level of background rejection 0, @] and/or 
to provide an unambiguous positive WIMP signal. The second strategy can 
be applied by searching for a correlation of the WIMP signal either with the 
motion of the Earth around the Sun, observed as an annual modulation of 
the number of events jsj, or with the direction of the Solar motion around 
the Galactic center (if, which happens to be roughly in the direction of the 
Cygnus constellation. The latter is generally referred to as directional detec- 
tion of Dark Matter and several project of detectors are being developed for 

this goal 53, Q E3, H Q. 

Directional detection of Dark Matter has first been proposed in [9| , highlight- 
ing the fact that a strong forward/backward asymmetry is expected in the 
case of an isothermal spherical Galactic halo, which could, in principle, be 
observed even with a poor angular resolution. Since then, several phe nomeno- 
logical studies have been performed 15^ W^L^ 18, 19, W 
Using an unbinned likelihood method 
tical tests on unbinned data 



18, 19, 20 
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161 . 171 or non-parametric statis- 



21, 22, 23, 24 



it has been shown that a few 
number of events O(10) is required to reject the isotropy, which is character- 
istic of the background event distribution. However, these methods do not 
allow to identify and characterize a WIMP signal and hence do not take full 
advantage of data from upcoming directional detectors [14 . 

In this Letter, we use a map-based likelihood analysis in order to extract 
from an observed recoil map the main incoming direction of the events and 
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its significance. In this way, the Galactic origin of the signal can be proved 
by showing its correlation with the direction of the Solar motion with respect 
to the Dark Matter halo. Hence, the goal of this new method is not to reject 
the background hypothesis, but rather to identify a genuine WIMP signal. 
This method is tested on realistic simulated cases defined by a low number 
of WIMP events O(10) hidden by a non-negligible background and measured 
with a rather low angular resolution ~ 15° (FWHM) [25|. This blind analysis 
is intended to be applied to directional data of any detector. 
The outline of this Letter is the following : after introducing the directional 
detection framework and phenomenology, the map-based likelihood analysis 
is presented and applied to a realistic simulated recoil map. Main results 
are presented and constraints are deduced in the (a n ,m x ) plane. Eventually, 
systematic studies are presented in order to show that this analysis tool gives 
satisfactory results on a large range of exposures and background levels. 



2. Directional detection framework 



2.1. Directional detectors 



Several directional detectors are being developed and/or operated : DRIFT 

|. A detailed overview of the 
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lOj, NEWAGE (111, MIMAC |12j, DM-TPC 
status of experimental efforts devoted to directional Dark Matter detection 
is presented in Directional detection of Dark Matter requires track re- 
construction of recoiling nuclei down to a few keV. This can be achieved with 
low pressure gaseous detectors |28] and several gases have been suggested : 
CF 4 , 3 He + C4H10 or CS2. Both the energy and the track of the recoiling 
nucleus need to be measured precisely. Ideally, recoiling tracks should be 3D 
reconstructed as the required exposure is decreased by an order of magnitude 
between 2D read-out and 3D read-out [20j. Sense recognition of the recoil 



track (head-tail) is also a key issue for directional detection 29|, |30 



2.2. Theoretical framework 

As a first application of the method, we focus on the simplest model for 
the Milky Way halo : the isotropic isothermal sphere, in which the WIMP 
velocity follows a Maxwellian distribution defined in the laboratory rest frame 
as 

(v + v G ) 2 ' 



1 



(27ra2) 3 / 2 



exp 



2af 



(1) 
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with a dispersion a v = v / v2 where t> = 220 km.s" 1 is the circular speed at 
Solar radius. As an illustration of the method, we consider a detector velocity 
equal to the tangential component of the Sun motion around the Galactic 
center v & = 220 ± 20 km.s -1 , neglecting the Sun peculiar velocity and the 
Earth orbital velocity about the Sun0. As a matter of fact, their contribution 



to detector velocity is less than the uncertainty on v & 18]. Using the Galactic 
coordinates (£, b), the WIMP velocity is written in the Galactic rest frame 
as : 

v = v (cos £ cos b x + sin I cos b y + sin b z) 

where x points towards the Galactic center, y in the direction of the Solar 
motion and z towards the Galactic north pole. The recoil distribution is then 
computed by generating random incident WIMP velocities from vf(v) and 
assuming an isotropic elastic scattering in the center of mass frame. Then, 
the recoil energy (En) in the laboratory rest frame is given by 

2v 2 m?m N „ , 

E * = 7 cos d * 2 

(m x + m N y 

with m x the WIMP mass, the mass of the target and 9r the recoil 
angle in the laboratory frame. Within this framework, the WIMP signal 
is expected to come from the Solar direction (£ Q = 90°, b Q = 0°) to which 
points the y axi^l. This happens to be roughly in the direction of the Cygnus 
constellation. 

In the following of the Letter, we will consider a form factor F(Er) taken 
to be equal to one. Indeed, using the Born approximation, in the case of a 
spin-dependent interaction, the form factor is given by the Fourier transform 



of a thin shell |26| leading to: F 2 (5 keV) = 0.99 and F 2 (50 keV) = 0.9 in 
the case of a 19 F target, justifying our approximation. Note that in the case 
of heavier targets, this approximation would be no longer valid. 

2.3. Recoil maps 

Recoil distributions are presented in Galactic coordinate maps, using the 
HealPix [27} tool in order to respect the spherical topology and to have an 



Obviously, when analysing real data, these two components of the detector velocity 
have to be considered in order to have an accurate analysis. 

2 However, when considering the peculiar velocity of the Sun, which is in the Galactic 
coordinate v Qp = ( 10.0, 7.3, 5.2) km.s -1 , the Solar motion direction is (£q — 87.5°, 6 Q = 
1.3°) 
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overall sky map with equal area bins. Figure [T] presents on the upper panel 
the theoretical WIMP flux and on the middle one, the theoretical recoil 
distribution evaluated with 10 8 WIMP-induced events, in the case of a 19 F 
target, a 100 GeV.c -2 WIMP and considering recoil energies in the range 5 
keV < Er < 50 keV. The lower bound of the energy range is due to the 
threshold ionization energy taking into account the quenching factor. As 
most of the WIMP events are concentrated at low recoil energy, an upper 
bound is chosen to limit the background contamination of the data. Indeed, 
following this framework, 70% of the recoils are between 5 keV and 50 keV 
and only 10% above 50 keV. Thus, increasing the upper bound would lead to 
a potentially weaker signal to noise ratio. It can be noticed that after scat- 
tering, the WIMP-induced recoil distribution, although obviously broadened, 
is still pointing in the Solar motion direction (£ Q ,b Q ). Thus, the expected 
WIMP signal is clearly anisotropic and directional detection should be able 
to distinguish a genuine WIMP signal from an isotropic background. 

2.4- WIMP signal characteristics 

In order to evaluate the evolution of the angular distribution of WIMP 
events as a function of the WIMP mass (m x ) and the recoil energy range 
[En), we first look at the ID angular spectrum, as it is a convenient repre- 
sentation of the recoil map. The latter is defined as the normalized fraction of 
events per solid angle as a function of the opening angle 7, corresponding to 
the angle between the y axis and the recoil direction. The result is twofold : 

• The upper panel of figure |2] presents the ID angular spectrum of WIMP 
events for three different values of m x : 10, 100, 1000 GeV.c -2 in the case 
of a 19 F target for a recoil energy in the range 5 keV < Er < 50 keV. We 
can notice that the WIMP event angular distribution clearly depends 
on the WIMP mass. Lighter is the WIMP, stronger is the angular 
anisotropy. Indeed, due to the finite energy range and the fact that low 
WIMP mass induce an energy distribution shifted to low energy, events 
above threshold are the one with the most directional feature (eq. (J2J ) • 
The directionality evolves quickly at low masses and very slowly for 
masses heavier than a hundreds of GeV.c -2 . Even if the WIMP is heavy 
(~ 1 TeV.c~ 2 ), the signal is still directional and then different enough 
from the background to be identified. Hence, directional detection 
presents the potential to distinguish a WIMP signal from background, 
for any WIMP mass. 
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• The lower panel of figure [2] presents, in the case of a 19 F target, three 
different ID angular spectra corresponding to three different recoil en- 
ergy ranges: 10 keV < E R < 20 keV, 30 keV < E R < 40 keV and 50 
keV < E R < 60 keV. We can notice that the whole angular distribution 
also depends on the recoil energy range and that greater is the recoil 
energy, stronger is the angular signature, see eq. §2§. 

The recoil map does also depend on the target mass rriN and as several 
directional detectors are being developed with different targets, we study 
the influence of these targets. Although their detection characteristics may 
be different (e.g. track length, drift velocity and straggling^]), for low mass 
targets (H, 19 F, 3 He) and at sufficiently low recoil energy when the form factor 
can be approximated to unity, equivalent directional signal can be found by 
adjusting the energy range for each target. Indeed, the directionality of 
the signal is encoded only in the cos 2 6 R term of eq. ((2]), then the angular 
distribution for a target of mass m Nl at a recoil energy E Rl is equivalent to 
the one of a target at Er 2 , using : 



Of course, this equivalence relation is no more valid for heavier targets when 
the form factor is far from unity. Hereafter, we consider a 19 F target with 
recoil energy in the range 5 < E R < 50 keV and m x = 100 GeV.c -2 (fig. 
[T]), but same WIMP-induced recoil distribution would be obtained for 3 He 
target with 1 < E R < 10 keV. 

3. Map-based likelihood analysis 

3.1. A realistic simulated measurement 

Lower panel of figured] presents a typical recoil distribution observed by a 
directional detector : 100 WIMP-induced events and 100 background events 
generated isotropically. These events are meant to be after data rejection 
based e.g. on track length and energy selection [l2j]. For an elastic axial cross- 
section on nucleon a n = 1.5 x 10~ 3 pb and a 100 GeV.c -2 WIMP mass, this 



3 This is the angular deviation due to scatterings of the recoiling nucleus with other 
nuclei from the gas. 



m N2 ( m x + m N] 



) 



2 



Er 2 = E Rl 



(3) 
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corresponds to an exposure of ~ 7 x 10 3 kg. day in 3 He and ~ 1.6 x 10 3 kg. day 
in CF 4 , on their equivalent energy ranges discussed in sec. 12.41 Low resolu- 
tion maps are used in this case (iV P i xe i s = 768) which is sufficient for a rather 
low angular resolution, ~ 15° (FWHM), expected for this type of detector 



and justified for instance by the straggling of the recoiling nucleus 25]. 3D 
read-out and sense recognition are considered for this study. 



3.2. Likelihood definition 

At first sight, it seems difficult to conclude from the recoil map of fig. [T] 
that it does contain a fraction of WIMP events pointing towards the direction 
of the Solar motion. A likelihood analysis is developed in order to retrieve 
from a recoil map : the main direction of the incoming events in Galactic 
coordinates (£, b) and the number of WIMP events contained in the map. 
The likelihood value is estimated using a binned map of the overall sky with 
Poissonian statistics, as follows : 

-^pixels 

Sf(m x ,\,£,b)= n P{[{l-X)B i + XS i {m x ;£,b))\M i ) (4) 
i=i 

where B is the background spatial distribution taken as isotropic, 5* is the 
WIMP-induced recoil distribution and M is the measurement. This is a four 
parameter likelihood analysis with m x , A = S/ (B+S) the WIMP fraction (re- 
lated to the background rejection power of the detector) and the coordinates 
(£, b) referring to the maximum of the WIMP event angular distribution. 
Hence, S(m x ; £, b) corresponds to a rotation of the S(m x ) distribution by the 
angles {£' = £ - £ Q , b' = b-b Q ). 

A scan of the four parameters with flat priors, allows to evaluate the 
likelihood between the measurement (fig. [1] bottom) and the theoretical dis- 
tribution made of a superposition of an isotropic background and a pure 
WIMP signal (fig. [T] middle). By scanning on £ and b values, we ensure that 
there is no prior on the direction of the center of the WIMP-induced recoil 
distribution. In order to respect the spherical topology, a careful rotation of 
the S distribution on the whole sphere must be done as follows. Given T^j 
the vector pointing on a bin Si, the following rotation is considered : 

f\ = R u {b')R s {£')f i 
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with u = Rz(£') x = u x x + u y y and Rg(b') the Rodrigues rotation matrix 
around an arbitrary vector u given by : 



where the angles are defined as follows : £' = £ — £ & , b' = b — b & . 

The events contained in the observed recoil map can be either WIMP- 
induced recoils or background events. A realistic recoil map from upcoming 
directional detectors cannot be background free. With this method, both 
components (background and signal) are taken into account and no assump- 
tion on the origin of each event is needed. Indeed, the observed map is 
considered as a superposition of the background and WIMP signal distribu- 
tions, and the likelihood method allows to recover A, the signal to noise ratio. 
The advantage is twofold : 

• First, background-induced bias is avoided. This would not be the case 
with a method trying to evaluate a likelihood on a map containing 
a fairly large number of background events considering only a pure 
WIMP reference distribution. 

• Second, the value of A allows to access the number of genuine WIMP 
events and consequently to constrain the WIMP-nucleon cross-section 
as presented in sec. 13.41 

It is worth noticing that the likelihood is performed on the whole angular 
distribution in order to maximize the information contained in the observed 
recoil map leading to more restrictive constraints on the four parameters. 
For instance, working on the forward/backward asymmetry is not sufficient 
to constrain the (£, b) parameters which are indeed the ultimate proof of the 
Galactic origin of the signal. 

3. 3. Results from a realistic recoil map 

The four parameter likelihood analysis has been computed on the sim- 
ulated map (fig. [1] bottom) and the full result is presented on figure [3j 
Marginalised distribution (diagonal) and 2D distribution (off-diagonal) plots 
of the four parameters m x , A, £, b are presented. The conclusion of the anal- 
ysis is twofold: 




u x u y {l — cos b') 
ib' + uj(l - cos b') 
u x sin b' 
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• Firstly, in the lower half part of figure [3] are presented the marginalised 
distributions of the parameters £, b and the 2D representation of Jzf in 
the (£, b) plane. These two parameters are well constrained and the first 
result of this map-based likelihood method is that the recovered main 
recoil direction is pointing towards (£ = 95° ± 10°, b = —6° ± 10°) at 
68 % CL, corresponding to a non-ambiguous detection of particles from 
the Galactic halo. This is indeed the discovery proof of this detection 
strategy. 

• Secondly, on the upper half part of figure |3] are presented the marginalised 
distributions of the parameters m x , A and the 2D representation of Jzf 
in the {m x , A) plane. The method allows to constrain A the WIMP frac- 
tion in the observed recoil map but not the WIMP mass. Constraining 
the WIMP mass is not the main point of this analysis tool. In fact, m x 
is set as a free parameter in order to show that the analysis is particle 
physics model independent. Indeed, the unknown value of m x does not 
affect the values of A, i and b due to the absence of correlations. Then, 
we can see on figure |3] that only m x < 10 GeV.c -2 is excluded. We 
can estimate the number of WIMP events as A^ wimp = A x N tot where 
Ntot = S + B follows a Poissonian statistic, and AA^ wimp is given by 
AiV wimp « A A x 

As a conclusion of this likelihood analysis, we can see that this method allows 
to determine from a realistic simulated recoil map (fig. [TJbottom) that it does 
contain a signal pointing towards the Cygnus constellation within 10°, with 
-Wwimp = 106 ± 17 (68%CL), corresponding to a high significance detection 
of Galactic Dark Matter. 

3.4- Constraining the elastic scattering cross-section 

A constraint in the (cr n , m x ) plane is then deduced from the marginalised 
=£?(A) distribution evaluated for each WIMP mass above 10 GeV.c" 2 . Using 
the standard expression of the event rate with a form factor F 2 (Er) taken 
equal to one and the local halo density p — 0.3 GeV.c -2 . cm -3 , the la and 
3cr CL contours are calculated. Figure @] presents the spin dependent cross- 
section on proton (pb) as a function of the WIMP mass (GeV/c 2 ). Results 
are presented in the case of a pure-proton approximation l31J and the proton 
spin content of 19 F has been chosen as < S p >= 0.441 [32|. The la and 3a 
contours deduced from the analysis of the simulated recoil map are presented 



9 



as shaded areas. It should be highlighted that these contours represent the 
allowed regions, as directional detection aims at identifying WIMP signal 
rather that rejecting the background. Exclusion limits from direct detection 

KIMS H, NAIAD 0, Picasso [38 



searches are presented 
and Xenon 10 
scope 



COUPP 35 
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Exclusion limit obtained with the Super-K neutrino tele- 
is also displayed. The theoretical region, obtained within the 
framework of the constrained minimal supersymmetric model, is taken from 
33j . Constraints from collider data and relic abundance Q x h 2 , as measured 
with WMAP 5-year data [34], are accounted for. The input value for the 
simulation is presented star. 

Such a result could be obtained, with a background rate of ~ 0.07 kg _1 day _1 
and a 10 kg CF 4 detector during ~ 5 months, noticing that the detector 
should allow 3D track reconstruction, with sense recognition down to 5 keV. 



4. Discovery potential 

In order to explore the robustness of the method, a systematic study has 
been done with 10 4 experiments for various number of WIMP events (iV wim p) 
and several values of WIMP fraction in the observed map (A), ranging from 
0.1 to 1. For a given cross-section, these two parameters are related re- 
spectively with the exposure and the rejection power of the offline analysis 
preceding the likelihood method. 

Figure presents on the upper panel the directional signature, taken as the 
value of a 7 = y/o^al, the radius of the 68% CL contour of the marginalised 
J?? [I , b) distribution, as a function of A. It is related to the ability to recover 
the main signal direction and to sign its Galactic origin. It can first be no- 
ticed that the directional signature is of the order of 10° to 20° on a wide 
range of WIMP fractions. Even for low number of WIMPs and for a low 
WIMP fraction (meaning a poor rejection power), the directional signature 
remains clear. From this, we conclude that a directional evidence in favor 
of Galactic Dark Matter may be obtained with upcoming experiments even 
at low exposure (i.e. a low number of observed WIMPs) and with a non- 
negligible background contamination (even greater than 50%). 
However, a convincing proof of the detection of WIMPs would require a di- 
rectional signature with sufficient significance. We defined the significance 
of this identification strategy as A/o^, presented on figure |5] (lower panel) 
as a function of A = S/(S + B). As expected, the significance is increasing 
both with the number of WIMP events and with the WIMP fraction, but 
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we can notice that an evidence (3a) or a discovery (5a) of a Dark Matter 
signal would require either a larger number of WIMPs or a lower background 
contamination. 

Using this map-based likelihood method, a directional detector may provide, 
as a first step, a Galactic signature even with a low number of WIMPs. For 
instance, a signal pointing towards Cygnus within 20° can be obtained with 
as low as 25 WIMPs with a 50% background contamination. For an ax- 
ial cross-section on nucleon of a Q = 1.5 x 10~ 3 pb, this corresponds to an 
exposure of 400 kg.day in CF 4 . In a second step, with an exposure four 
times larger, the directional signature would be only slightly better (10°) but 
the significance would be much higher (~ 7a) and the detection much more 
convincing. 



5. Conclusion 

We have presented a powerful statistical analysis tool to extract infor- 
mation from a data sample of a directional detector in order to identify a 
Galactic WIMP signal. As a proof of principle, it has been tested within 
the framework of an isothermal spherical halo model. We have shown the 
feasibility to extract from an observed map the main incoming direction of 
the signal and its significance, thus proving its Galactic origin. Evidence in 
favor of Galactic Dark Matter may be within reach of upcoming directional 
detectors even at low exposure. In a second step, with increasing exposure 
and by using an extended method currently under development, Dark Matter 
properties (triaxility, halo rotation, WIMP mass) may also be constrained. 



For this first study of directional detection as a strategy to discover Galac- 
tic Dark Matter, the standard halo model has been considered. It may be 
noticed that recent results from high resolution numerical simulations seem 
to be in favour of non-smooth WIMP velocity distributions 4l|, |42] . However, 



as their resolution is many orders of magnitude larger than the scale of the 
ultra-local Dark Matter distribution probed by current and future detectors, 
this result is not fully relevant for direct detection. Even if applied at ultra- 
local scale, the result from 42] is that high velocity WIMPs (v > 500 km.s -1 ) 
could deviate from the Solar motion direction (£ Q , b Q ) within about 10° due 
to the presence of subhalos, tidal streams or an anisotropic velocity distri- 
bution. But, as the minimal speed is 130 km.s" 1 in this study and as the 
Maxwellian distribution of WIMP velocity is peaked around 300 km.s -1 the 
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contribution of high speed WIMPs to the WIMP flux is negligible. For in- 
stance, this high speed feature would be present in data from detector with 
a higher energy threshold. In such case, this will mildly affect the result pre- 
sented in this Letter by shifting the main recoil direction recovered from the 
likelihood analysis by a few degrees. However, the discovery proof will still be 
reached as long as the WIMP signal is still directional, thus remaining differ- 
ent from the background. Effect of non-standard halo models (triaxial, with 
stochastic features or streams) will be addressed in a dedicated forthcoming 
paper. 



Note added after submission 

While this work was being refereed, a related work from A. M. Green and 
B. Morgan 43| appeared on the arXiv. Their paper is complementary to this 
one as it is estimating the number of WIMP events required to confirm the 
direction of Solar motion (£ & ,b & ) as the median inverse recoil direction at 
95% CL. Then, both this work and the one from [43[ highlight the fact that 
directional detection is a key issue to clearly identify a positive Dark Matter 
detection. 



12 



References 

[1] V. C. Rubin, A. H. Waterman and J. D. P. Kenney, The Astronomical 
Journal 118 (1999) 236-260 

[2] D. Clowe et ai, Astrophys. J., 648 (2006) L109 

[3] E. Komatsu et ai, Astrophys. J. Suppl. 180 (2009) 330 

[4] G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267 (1996) 
195 

[5] G. Servant and T. M. P. Tait, Nucl. Phys. B 650 (2003) 391 
[6] Z. Ahmed et ai, arXiv:0912.3592l 

[7] for a recent review, see e.g. Proc. of Identification of Dark Matter 2008 
(IDM 2008), Stockholm, Sweden, Aug 2008 

[8] K. Freese, J. Frieman and A. Gould , Phys. Rev. D 37 (1988) 3388 

[9] D. N. Spergel, Phys. Rev. D 37 (1988) 1353. 

[10] G. J. Alner et a/., Nucl. Instr. Meth. A 555 (2005) 173. 

[11] K. Miuchi et ai, Phys. Lett. B 654 (2007) 58 

[12] D. Santos et ai, J. Phys. Conf. Ser. 65 (2007) 012012 

[13] G. Sciolla et al. larXiv:0806.2673l 

[14] S. Ahlen et ai, International Journal of Modern Physics A25 (2010) 1-51 

[15] C. J. Copi and L. M. Krauss,, Phys. Rev. D 63 (2001) 043507 

[16] C. J. Copi, L. M. Krauss, D. Simmons-Duffin and S. R. Stroiney, Phys. 
Rev. D 75 (2007) 023514, 

[17] C. J. Copi and L. M. Krauss, Phys. Lett. B 461 (1999) 43 

[18] B. Morgan, A. M. Green and N. J. C. Spooner, Phys. Rev. D 71 (2005) 
103507, 

[19] B. Morgan and A. M. Green, Phys. Rev. D 72 (2005) 123501 

13 



[20] A. M. Green and B. Morgan, Astropart. Phys. 27 (2007) 142, 

[21] A. M. Green and B. Morgan, Phys. Rev. D 77 (2008) 027303 

[22] O. Host and S. H. Hansen, JCAP 0706 (2007) 016 

[23] J. D. Vergados and A. Faessler,, Phys. Rev. D 75 (2007) 055007 

[24] M. S. Alenazi and P. Gondolo, Phys. Rev. D 77 (2008) 043532 

[25] D. Santos et al, in preparation 

[26] J. D. Lewin & P. F. Smith, Astropart. Physics 6 (1996) 87-112 

[27] K. M. Gorski et al, Astrophys. J. 622 (2005) 759 

[28] G. Sciolla and C. J. Martoff, New J. Phys. 11 (2009) 105018 

[29] D. Dujmic et al, Nucl. Instrum. Meth. A 584 (2008) 327 

[30] S. Burgos et al. larXiv:0809.1831l 

[31] D. R. Tovey et al, Phys. Lett. B 488 (2000) 17 

[32] A. F. Pacheco and D. Strottman, Phys. Rev. D 40 (1989) 2131-2133 

[33] R. Trotta, F. Feroz, M. P. Hobson, L. Roszkowski and R. Ruiz de Austri, 
JHEP 0812 (2008) 024 

[34] J. Dunkley et al, Astrophys. J. Suppl. 180 (2009) 306 

[35] E. Behnke et al, Science 319 (2008) 933 

[36] H. S. Lee et al, Phys. Rev. Lett. 99 (2007) 091301 

[37] G. J. Alner et al, Phys. Lett. B 616 (2005) 17 

[38] S. Archambault et al, Phys. Lett. B 682 (2009) 185 

[39] J. Angle et al, Phys. Rev. Lett. 101 (2008) 091301 

[40] S. Desai et al, Phys. Rev. D 70 (2004) 083523 

[41] M. Vogelsberger et al, Monthly Notices of the Royal Astronomical So- 
ciety 395 (2009) 797-811 

14 



[42] M. Kuhlen et al, JCAP 1002 (2010) 030 

[43] A. M. Green and B. Morgan, Phys. Rev. D 81 (2010) 061301 



15 




0.0 ^^^^^^ht ^^^^H 4.0 Number of events 



Figure 1: From top to bottom : WIMP flux in the case of an isothermal spherical halo, 
WIMP-induced recoil distribution and a typical simulated measurement : 100 WIMP- 
induced recoils and 100 background events with a low angular resolution. Recoils maps 
are produced for a 19 F target, a 100 GeV.c -2 WIMP and considering recoil energies in 
the range 5 keV < Er < 50 keV. Maps are Mollweide equal area projections. 
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Figure 2: Normalized angular spectra of the WIMP signal, i.e. fraction of events per solid 
angle as a function of the opening angle 7, defined as the angle between the y axis and 
the recoil direction. All plots are in the case of a 19 F target. Upper panel presents the 
normalized angular spectra for three different values of m x and a recoil energy range 5 keV 
< Er < 50 keV. Lower panel presents the normalized angular spectra for three different 
recoil energy ranges and with m x = 100 GeV.c -2 . Black solid line presents the isotropic 
background spectrum. 
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Likelihood analysis 



WIMP mass: m_ t > 10 GeV.c" 2 
WIMP fraction: X = 0.53 ± 0.085 (1oCL) 
Galactic latitude: I = 95° ± 10° (1a CL) 
Galactic Longitude: b = -6°± 10° (1 a CL) 
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Figure 3: Marginalised distribution (diagonal) and 2D distribution (off-diagonal) plots of 
the four parameters m x , A, £, b from the analysis of the simulated map (fig. [T] lower panel). 
Blue squares represent the calculated values of and orange triangles refer to the input 



value of each parameter: 



100 GeV.c -2 , A 



0.5,^ = 90°, 6 = 0°. 
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Figure 4: Spin dependent cross-section on proton (pb) as a function of the WIMP mass 
(GeV/c 2 ). Results are presented in the case of pure-proton approximation (3l|. The theo- 
retical region, obtained within the framework of the constrained minimal supersymmctric 
model, is taken from [33(. Constraints from collider data and relic abundance are ac- 



counted for. Exclusion limits from direct detection searches are presented : COUPP [35j], 
KIMS NAIAD Picasso [H| and XenonlO j3§|. Exclusion limit obtained with the 
Super-K neutrino telescope [Io| is also displayed. The allowed regions obtained with the 
example map shown are presented with shaded area, with la and 3u CL contours. Input 
value for the simulation is shown with a star. 



19 




0.8 1 

X (wimp fraction) 



15 



10- 



- N wimp = 25 

N vvimp = 50 

* N wimp = 100 



0.2 



-k' 
Y - 



T" 



0.4 



0.6 



X (wimp fraction) 



Figure 5: Upper panel presents the directional signature cr 7 (in degrees), the radius of the 
68 % CL contour of the Jzf (f, b) distribution, as a function of the WIMP fraction. Lower 
panel presents the significance (A/<ta) as a function of the WIMP fraction A = S/(S + B). 
Results are produced for a 19 F target, a 100 GeV.c" 2 WIMP and considering recoil energies 
in the range 5 keV < E R < 50 keV. 
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